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Abstract 


An important consideration in the use of manipulators in Microgravity environments is the 
minimization of the base reactions, i.e. the magnitude of the force and the moment exerted by the 
manipulator on its base as it performs its tasks. One approach which has been proposed and 
implemented is to use the redundant degrees of freedom in a kinematically redundant 
manipulator to plan manipulator trajectories to minimize base reactions. In this paper we 
develop a global approach for minimizing the magnitude of the base reactions for kinematically 
redundant manipulators which integrates the Partitioned Jacobian method of redundancy 
resolution, a 4-3-4 joint-trajectory representation and the minimization of a cost function which 
is the time-integral of the magnitude of the base reactions. We also compare the global approach 
with a local approach developed earlier for the case of point-to-point motion of a three degree- 
of -freedom planar manipulator with one redundant degree-of-freedom. The results of study 
show that the global approach is more effective in reducing and smoothing the base force while 
the local approach is superior in reducing the base moment. 
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1. Introduction 

Recently there has been a lot of interest in the use of kinematically redundant manipulators for various 
applications like obstacle avoidance and singularity avoidance [2-4]. In this paper we want to utilize the well-known 
idea that kinematic redundancy provides the analyst or designer with additional degrees of freedom which can be 
exploited for "useful" purposes. Specifically we are interested in the planning of trajectories for manipulators which 
are used in space applications [1] and typically operate in microgravity environments. 

An issue of considerable importance for manipulators used in Microgravity environments is the minimization 
of the magnitudes of the dynamic reaction force and moment exerted by the manipulator on its base as it performs 
its task. One reason for minimizing, and if possible eliminating base reactions, is that the base force and base 
moment could disturb other tasks or experiments in the vicinity. 

In this paper we address the issue of using kinematic redundancy to plan joint-space trajectories which 
minimize the reactions transmitted by the manipulator to the base. A local approach to this problem is described and 
implemented in [8-9]: in this approach a performance index or cost function was minimized at each time segment. 
One of the problems with the local approach was that it led to undesirable "peaks" in the base reactions as shown in 
Fig.(8). In this paper, we propose and implement a global approach which is based on minimizing an integral 
performance index over the entire end-effector trajectory. 

The rationale which underlies our approach is as follows. The trajectory-planning problem for a manipulator 
reduces to the problem of solving the inverse-kinematic problem for the joint-variables (or "joint-space" trajectory ) 
given the specified trajectory in the task space. In the case of a (kinematically) redundant manipulator the inverse- 
kinematic problem has an infinite number of solutions. The basic idea is to pose an optimization problem of 
minimizing the base reactions in order to obtain joint-space trajectories (which minimize base reactions). In this 
paper we develop a global approach which integrates the Partitioned Jacobian method of redundant resolution, a 
4-3-4 joint-trajectory representation (descrilred in Section 3.3) and the minimization of a cost function which is the 
time-integral of the magnitude of the base reactions. We also compare the global approach with a local approach 
developed earlier for the case of point-to-point motion of a three degree-of-freedom planar manipulator with one 
redundant degree-of-freedom. 

The contents of the paper are organized as follows. In Section 2 we develop a statement of the problem to be 
addressed in this paper. In Section 3, we show how tire Partitioned Jacobian method of resolving kinematic 
redundancy can be used to develop a Global approach for minimizing the magnitude of the base reactions. In this 
section we also briefly describe the local approach used in [8,9] since we are interested in comparing the two 
approaches. The global and local approaches are then applied (Section 4) to obtain the optimal joint trajectories 
which minimize base reactions for a redundant three degrees of freedom planar manipulator (with 1 redundant 
degree of freedom for its intended task of moving from point-to-point in the workspace). A simple sensitivity study 
is also carried out in order to get some feel for the "goodness" of the optimized trajectory. The results of the case 
study (of Section 4) are then discussed in some detail in Section S. Finally, in Section 6 we summarize the work and 
draw some conclusions. 

2. Statement of the Problem 

Consider an m degrees of freedom manipulator with m revolute or prismatic joints. Let q £ denote the joint 
variable at joint i, the joint variable of each joint is either an angle 6; (revolute joint) or a distance dj (prismatic 
joint). For an m degrees of freedom manipulator, the joint displacement vector q can be defined as 

q - [ qj. q 2 . (1) 

Let E be an arbitrary point on the end-effector as shown in Figure 1. The position of E can be represented with 
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reference to the coordinate frame {U} fixed in the base by 

r = [ x y z ] T . (2) 

The kinematic equations relating the end-effector position r and the joint variable vector q can be expressed as 
follows: 

r(0 = g(q(0). (3) 

where g is a vector of functions which are nonlinear in q(t). 

The trajectory planning problem for manipulators is to determine the joint variables q(t) for a specified 
end-effector trajectory r(t). For a redundant manipulator, the number of degrees of freedom, m is more than the 
minimum number of degrees of freedom required to perform a task which is denoted by n. Therefore, from equation 
(3) we can see that there are an infinite number of joint space solutions q(t) for a specified end-effector position, 
r(t). This means that we can take advantage of the redundancy of (kinematically) redundant manipulators to 
optimize some suitable dynamic criterion. The optimization of this dynamic criterion will then, in general, yield a 
unique (joint space) solution for the joint variables q(t). 

An important problem which arises in the use of manipulators in space is the minimization of the magnitudes 
of the reaction force and moment transmitted to the base of the manipulator as the end-effector of the manipulator 
follows a specified trajectory. Let F Q and N c , respectively, denote the force vector and moment vector (see Figure 
1) transmitted to the base of the manipulator and define 

F b T =t F o TN o T l- W 

We will call F b the base reaction vector. The dynamic equations relating the motion of the manipulator to the force 
and moment transmitted to the base (in a microgravity environment) can be represented in the following form: 

F b = M b (q)q+V(q,q), (5) 

where M b (q) is a mass matrix and V(-) denotes all terms which are nonlinear in q. (In Appendix 2 we give 
expressions for F b , M b and V for a planar 3 d.o.f. manipulator). 

Since we are interested in minimizing the base reactions, we must select a cost function B which is an 
appropriate function of the base reaction vector F b , i.e. 

B=K F b ). (6) 

Since F b is a function of (q, q, q), B is also a function of (q, q, q). Therefore minimization of the magnitude 
of the base forces and base moments reduces to determining the joint space trajectory q(t) which will mimimize the 
cost function B. 

We are now in a position to state the objectives of this paper 

1. Given a specified trajectory of the end-effector, use kinematic redundancy to develop an approach for 
planning the corresponding joint-space trajectories q(t) in order to mimimize a cost function which is a 
measure of the time-integral of the magnitude of the base reactions. This approach is called the Global 
Approach. 

2. Compare the global approach with the local approach developed in [8,9]. 

3. Study the sensitivity of the (global) cost function with respect to (small) changes in the joint 
trajectories. 

In the next section we develop the global approach and very briefly describe the local approach. 
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3. Description of Approach 

The general problem of obtaining unique inverse kinematic solutions for redundant manipulators is termed 
"redundancy resolution". Even though there are a variety of approaches (which are mostly based on the pseudo- 
inverse of the manipulator Jacobian) proposed for resolving the redundancy of redundant manipulators, relatively 
few approaches [5,6,7,8,9] address the problem of utilizing redundancy for optimizing dynamic criteria such as 
torque optimization [5,6,7] or base reaction minimization [8,9]. Hollerbach and Suh [5] investigated a few local 
approaches based on the pseudo-inverse for optimizing joint torques. The results of their study show that the local 
approach leads to "unstable" solutions for some trajectories. Nakamura [7] proposed a pseudo-inverse approach 
based on the Pontrayagin’s Maximum Principle for globally optimizing the energy used in performing task. Suh and 
Hollerbach [6] investigated the problem of torque optimization using local and global approaches based on the 
pseudo-inverse. It was found that the global approach does not have a stability problem and outperforms the local 
approach [5] at all times. However, the global approach proposed in [6,7] is complex to formulate and solve, and 
limits on joint angle, joint velocity, and torques cannot be included in the optimization problem. The shortcomings 
of the pseudo-inverse approach in general are clearly given in [14]. 

In contrast to the pseudo-inverse approach, we describe an alternative approach, based on partitioning the 
manipulator Jacobian matrix [4], which is particularly well suited to the base reaction minimization problem. The 
first application of the Partitioned Jacobian approach to the base reaction minimization problem was a local 
approach described in [8,9]. The primary aim of this paper is to develop and implement a global approach which 
utilizes the Partitioned Jacobian method of resolving redundancy. The global approach is described in Section 3.3. 
Since we also want to compare the global approach with the local approach, the local approach is described briefly 
in Section 3.2. The Partitioned Jacobian method of resolving redundancy is explained in Section 3.1. 

3.1. Redundancy Resolution 

Consider an m degree-of-freedom redundant manipulator. If the velocity v of a reference point on the 
end-effector is represented by an (nxl) column vector, then we can write the following relationship between v and 

the (mxl) column vector q of rates of change of joint variables: 

v = J(q)q, (7) 

where J is an (nxm) manipulator Jacobian matrix, (q will be referred to as the joint "velocity" vector for short) 

Differentiating Eq. (7) with respect to time, the acceleration of the end-effector is 

v = j q + J q. (8) 

For redundant manipulators, the Jacobian matrix, J is not a square matrix. Therefore for a given end-effector 

velocity v, the Jacobian cannot be inverted to obtain the corresponding joint "velocity" vector q. In order to obtain 
joint space solution, we use the Partitioned Jacobian approach developed in [4]. Using this approach, the joint 
variable vector can be partitioned into, 

q r - [ q„ T q r T ). (9) 

where q r an (m-n)xl column vector, is called the redundant joint vector and q nr , an nxl column vector, is called the 
nonredundant joint vector. 


Eq. (7) can then be written in the form 

V = Jn, + *r <lr> 


( 10 ) 
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where J nr is an nxn matrix corresponding to and J r is an nx(m-n) matrix corresponding to q r . 

We have therefore partitioned the Jacobian J as follows 

J = [ J„r Jr 1 (H > 

Commenr .lt is always possible to find a vector of n non-redundant joint variables q^ such that J nr is an 
invertible (i.e. non- singular) matrix except at a finite number of singularities. 

Eq. (8) can be written as follows: 

V - K + 1 <ir + V + Jr V < 12 ) 

Rearranging Eqs. (10) and (12), we have 

<u, - J„-> ( * - J, q, ) (13) 

* J.,"' ( V - j, q, - j„ q„ - J, q r ). (14) 

In the approaches described below, one specifies the redundant joint variable vector q r (t) to be a convenient 
(vector) function of time from which q r (t) is readily obtained. Equation (13) can then be used to determine the 
non-redundant joint variable velocity vector q nr (t); from the latter q nr (t) can be obtained by integration. 

3.2. Local Approach 

In the local approach, the total time of the task is divided into a number of time segments. For convenience, 
assume that the duration for each time segment is the same. If q ro is the joint velocity vector of the redundant joints 

at the start of any time-segment At and q r is the joint velocity vector of the redundant joints at the end of the time 
segment, then we can write 

q r = q„+Aq r , (15) 

where Aq r can be approximated by a linear combination of shape-function f(t), i=l,2, Jc as follows 

A V « c, (r) + c l2 /*«) + ... + c 1(k /*(r) 

A q r 2 - c 21 / (0 + + ... + c^f) 

+ + + (16) 

where Aq t \ i=(l^ (m-n)), is the 1 th element of Aq r The use of shape functions in optimization of torques of 

non-redundant manipulators is described in [13]. 

If C denotes the matrix of parameters Cy, i=(l,..,(m-n)), j=(l^,..Jc), and f(t) is a column vector whose i lh 
element is f^t), then Eq.(16) can be written in a more compact form 


Aq r = Cf(t). 

Combining Eq.(15) and Eq.(17) we obtain 


( 17 ) 
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q r = q ro +Cm < 18 ) 

For convenience, f*(t) is chosen to be polynomials t 1 ' 1 , (i=l,2,..,k) and therefore f(t) is given by 

f(t) = [ 1 t t 2 ... t k - 1 ] T , (19) 

where k is the degree of the polynomial selected for q r . 

Substituting Eq.(18) in Eqs.(13-14), they become 

q„, = V 1 [ v - J t (q„+C0] (20) 

q„ - J„ 1 (vj/q^CO-U., 1 [v-J r (q ro .Cf)].J,Cf ) . (21) 


Now, we can compute the values of q, q, and q as functions of the elements c- of the C matrix for each time 
segment. For an arbitrary manipulator configuration, the base reactions (F 0 and Nj) at the end of each time step are 
therefore functions of the elements of the C matrix. In Section 3.4 (see Eq. 37) we will present a local cost function, 
B(C) that is a measure of the magnitudes of the base force and base moment at the end of each time segment. Using 
that cost function we can pose the inverse kinematic problem as an minimization problem which yields the optimal 
C matrix for each time segment. For further details refer to [8,9]. 


3 J. Global Approach 

In contrast to the local approach where the magnitudes of the base reactions at the end of each time segment is 
optimized, the global approach described below optimizes a cost function which is a time-integral of the magnitudes 
of the base force and base moment over the whole time interval, x f . 

Using the redundancy resolution scheme developed in Section 3.1, the redundant joint trajectories are used in 
the optimization of appropriate dynamic criteria such as the minimization of the magnitudes of the base reactions. 
The redundant joint velocity vector representation, Eq. (18) used in the local approach would be difficult to 
implement for a time integral cost function. In the local approach, for each time segment we have (m-n)xk 
parameters (i.e. the elements c^ of the matrix C). If the total time is divided into / time segments, then we would 
have lx(m-n)xk parameters to determine for global approach. To keep the number of parameters required for the 
representation of the redundant joint variables small, we use the simpler representation described below. 

One of the intuitive requirements for the redundant joint trajectories is that the joint trajectories have to be 
smooth. Also, the angular velocity and acceleration at both of the initial and final points of the redundant joint 
trajectories should be zero. One way of representing this family of curves is the use of higher order polynomials 
(splines) to satisfy the velocity and accelerations boundary conditions. Since very high order polynomials exhibit 
unpredicted behaviors, a popular method is to break a joint trajectory into a number of segments which can be 
represented by lower order polynomials such as third and fourth order ones. These segments will satisfy the 
boundary conditions and are continuous from one segment to another. The representation considered in this paper is 
a three-segment method which consists of two fourth-order polynomials far the initial and final segments and a third 
order polynomial for the transition segment. This method is known as the 4-3-4 joint trajectory representation. The 
details of this representation can be found in [10]. As shown in Figure 2, a typical 4-3-4 joint trajectory for each 
joint consists of three segments:an initial section (0j(t)), a transition section (0 2 (t», and a final section (0 3 (t)). 
Basically, curves described by this method are functions of the initial and final joint points (0j and 0 f ), the time 
period of the transition segment^), the total time (t f ) and the angular displacement traversed during the transition 
period (62). The expressions for 0j(t), 0 2 t), and 0 3 (t) are given in Appendix 1. 
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3J.1. Three-parameter Joint Trajectory Representation 

The 4-3-4 joint trajectory can be simplified by considering the following normalized parameters: (a) time 
factor, p (0 < p < 1 ) 


P = 


Tj +T2 + T3 


where tj’s are time intervals for the three segments (see Fig. 2); 


( 22 ) 


(b) angular displacement factor, a 

% 

""WV 

where Sj’s are angular displacements traversed during the ith time interval (see Fig. 2), i=l,2,3. 


(23) 


For x j - x 3 and 8 j = 83 , we have the following relationships: 



(24) 

8 2 =0 / O « 

(25) 

83 -V* 

(26) 

Xj^p' 

(27) 

x 2 =X/P 

(28) 

ts-'VP'. 

(29) 


where 0 fo = 0 f - 0 it x f is the total time of the task, a' = (1 - a )/2 and P' = (1 - p)/ 2 . 

Note: the above variables, defined by Eqs. (24) through (29) are functions of a, x f , p,0 it and 0 f . In general, for 
the redundant manipulator inverse kinematic problem, the initial angular position^) for the redundant joint variable 
and the total time of the task (x f ) are usually specified. Therefore, there are only three unknown paremeters, namely 
a, 0 f , and p for each degree of redundancy. Using these parameters, we obtain a three-parameter joint trajectory 
representation (based on the 4-3-4 joint trajectory representation) for the redundant joint variables. The j 111 element 
of the redundant joint velocity vector, qj can be expressed as 

*/(0-W f° r >1.23 (30) 

where 0 k i(t) denotes the polynomial for the Jfc* segment, (k=l,2,3), of the /** redundant joint variable, j=l, 2 ,..(m-n) 
and t^s are shown in Fig. 2. 

Figure 3 illustrates the curves generated by different combinations of a and p. Note that these curves are 
smooth and satisfy the boundary conditions. 


In the next section, we show how the three-parameter joint trajectory representation for redundant joint 
variables can be used to find a joint space solution (i.e. inverse kinematic solution) for the redundant manipulators. 


3.3.2. Global Optimization Problem Formulation 

For a redundant manipulator with (m-n) degrees of redundancy, the application of the three-parameter scheme 
to represent each of the (m-n) redundant joint variables yields q r as a function of 3x(m-n) number of parameters, o’, 
pi, and 0^, j= 1 , 2 , ..(m-n), which need to be determined. Using Eq. (13), q nr can also be determined as a function of 
these 3x(m-n) parameters. Let a, §, and 0f denote (m-n)xl column vectors whose /** element is denoted, 
respectively, by oi, P>, and 0 fi. 
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The cost function B which is a function of the base reactions and therefore a function of q, q, q, can 
expressed as a function of a, jj, and Of (see procedure below for details). Therefore the optimization problem can 
posed as follows: 


minimize B (a, J3, Of) 
subject to the constraints 

«La<<V<«U 

0<o>< l,and 

0 < [Ji < l,j=l,2,..(m-n). (31) 

For the global approach B is a time-integral cost function defined in Section 3.4 (Eq. 35) below, which is a measure 
of the magnitude of the base reactions. 6> min and 0i mtx are selected based on the information on the joint limits. 

Fig. 4 depicts a simplified flow-chart which captures the essence of the the global approach. The procedure 
for the global approach, corresponding to the flowchart in Fig. 4, is as follows: 

1. Initialize the parameters vectors a, §, and Of . 

2. Using the 4-3-4 redundant joint representation (Eq. (30) and Appendix 1) obtain q r , q r , and q r as 
functions of a, J3, and Of for the time interval x f . 

3. Using Eq. (13) and (14), respectively, to determine q nr and q nr , as a function of a, J3, and Of. 

4. Obtain the non-redundant joint variable vector q nr by integrating q nr (obtained in Step 3). We now 
have q, q, and q as a function of a, J3, and Of. 

5. Use the dynamic equations (Eq. (5)) and the results of Step 4 to determine the base reaction vector F b , 
as a function of a, §, and Of. 

6. Define a suitable time-integral cost function B of the base reactions F b (see Sections 3.4 and 4.2). 

7. Combine the results of steps 5 and 6 to obtain the cost function B as a function of a, J3, and Of. 

8. Apply the Hooke and Jeeves optimization method [12], to update the parameter vectors a, J3, and Of. 

/\ A A 

9. Repeat Steps 2 through 8 until the parameters vectors converge to their optimal values a, ]3, 0 f . 

Once the optimal values of a, JJ, and Of have been obtained, the optimal joint space trajectories vector q(t) are 
readily determined. 


3.4. Cost Function 

The selection of a suitable cost function is crucial in obtaining a joint space solution that minimizes base 
reactions. In this section, we discuss the issue of selecting a proper cost function which can be used as a measure of 
the magnitudes of the base reactions (F 0 and N^. 

In general, the magnitude of the base force is given by 

F o = ( F oi 2 + F o2 2 + F o3 2 ) 1/2 - (32) 

and the magnitude of the base moment is 

N o -<N o1 2 + N* 2 + N 03 *)»«. (33) 

We introduce two weighting factors (wj and w 2 ) that enable us to place appropriate weights on the base force 


£ £ 
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Fig. 4 Basic Algorithm For the Global Approach 
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and base moment components. The cost function B can be expressed as 

B=w x F 2 +w 2 N 2 , for the local approach , and (34) 

B= f / h ' 1 F a 2 + w 2 N 2 dt, for the global approach. (35) 


0 

w 2 l . 

We can write Eqs. (34) and (35) in a more compact matrix form as follows: 

B = F b T QF b 
and 

B=f T/ F/ Q F b dt. 

*0 

The weights and w 2 must be carefully chosen to : (l)appropriately scale the base reactions with respect to 
a reference set of base reactions and (2)ensure that all quantities in the cost function are dimensionally 
homogeneous. The choice of appropriate weights will be dealt with in Section 4.2. 

4. Illustrative Example 

In this section we apply the basic approaches described in the paper to minimizing the base reactions of a 
three degree-of-freedom planar manipulator. The task is a point-to-point motion of the end-effector from a specified 
initial position to a specified final position. Using the above approaches, we would like to plan a trajectory in the 
joint space to minimize the magnitudes of the base reactions as the manipulator moves from the initial position to 
the final position. 

In this problem, for the planar manipulator, m=3; since we are interested in the position of the end-effector, 
n=2. Therefore we have (m-n)=l degree of redundancy. The redundant link of this problem is chosen to be the first 
link which is attached to the base (q r «q t ). 

The base reactions dynamic equations of the 3 d.o.f. planar manipulator are in Appendix 2 and the link 
dimensions and mass properties are given in Appendix 3. 

We study the following four cases: 

1. The minimization of the time-integral of the magnitude of the base moment only. 

2. The minimization of the time-integral of the magnitude of the base fore only. 

3. The minimization of the time-integral of an appropriately weighted combination of the magnitudes of 
the base force and base moment. We use the results of cases (1) and (2) to show how the weights Wj 
and w 2 should be chosen. 

4. A sensitivity analysis of the sensitivity of the cost function of Case (3) with respect to small changes in 
(the parameters of) the joint space trajectories. 


(36) 

(37) 

(38) 




Define 


Q= 


w , 1 
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4.1. End-effector Trajectory 

The distance measured from the initial position along a straight line which connects the initial position and the 
final position can be expressed as 

dQ^bit-^sinA, (39) 

ZK Ty 

where b is an arbitrary constant and x { is the total time of the task. The speed of the end effector is given by 

v(t)-b(l-cos^-)). (40) 

x f 

Fig. 5 shows the desired velocity and acceleration profiles of the end-effector. 

The initial position of the end effector is (x,y) = (0.0, 1.2071) (m) which corresponds to an initial manipulator 
configuration, q T =(45°, 45°,45°]. The desired final position is (x,y) = (0.3536, 0.8535) (m). The following values 
are chosen for the trajectory: 

x f = 2.0s 
b « 0.25 
d(t f ) « 0.5 m 
A ts 0.1 s. 


42. Cost function [Case Studies 1-4] 

To minimize the base reactions properly, it is important to select a suitable weighting matrix Q for the cost 
function. Therefore, the choices of the two weighting factors (wj and w 2 ) are crucial. To understand the effects of 
these weighting factors on the base reactions, three cases of different Q matrices are studied. We only show the cost 
functions for the global approach. The cost function for the local approach can be obtained by Eq. (37). 

(a )Case Study 1 


In this case study we select a cost function that only minimizes the magnitude of the base moment in order to 
determine' how small the magnitude of the base moment can be made. We simply choose wj= 0 and w 2 -l in Eq. 
(36) to obtain the weighting matrix 

" 0 <f 


Qi= 


i 


and the corresponding cost function (from Eq. 38) 


(41) 


(42) 


(b) Case Study 2 


In this case, we want to know the extent to which the magnitude of the base force component can be 
minimized by choosing w t =1 and w 2 » 0 in Eq. (36) to obtain the matrix 

fl ol 


q 2 = 

o o 

and the corresponding cost function (from Eq. 38) takes the form 


(43) 



( 44 ) 
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(c) Case Study 3 

In general, we want to minimize the magnitudes of both the base force, F 0 and the base moment, N 0 
transmitted to the base. Since F 0 and N 0 are different physical quantities, we require a weighting matrix Q to 
non-dimensionalize and to appropriately scale the base force and base moment components. The results of Case 
studies 1 and 2 are important in this regard. 


The average value of the base force (F^ transmitted to the base in Case Study 1, F Jvg is a measure of how 
small we can make the base force. Similarly, the average value of the moment transmitted to the base in Case Study 
2, N Jvg is a measure of how small we can make the magnitude of the base moment. We therefore choose, 


-ttfr * >■ 

u Javg 1 avg 


( 45 ) 


The weighting matrix therefore is given by 

r ^ 


where 


Qs= 


wjl 0 
0 W 2 I , 


( 46 ) 


1 

W.- r 

1 F 2 

avg 

W,= r. 

N 2 

avg 


( 47 ) 

( 48 ) 


Note that the above choices of wj and w 2 simultaneously achieve the desired scaling and non 
dimensionalization. 


(d)Case Study 4 ( Sensitivity ) 

In practice, it is impossible to realize the exact optimal joint trajectories due to the errors in the control system 
and noise and uncertainties in the system. Therefore, it is important to examine how sensitive the base reactions are 
to the deviations of the joint trajectories from the optimal solution. The global approach provides us a simple means 
of studying this issue. One way of studying this issue is to introduce perturbations about the optimal values of a, p, 
and 6f obtained from Case Study 3. The following three cases in each of which one of the three parameters is 
increased by 20% were studied: 

case a : the optimal value of a, viz & is replaced by 1.2 a 

case b : the optimal value of p, viz fi is replaced by 1.2 0 

case c :the optimal value of 0 f , viz § f is replaced by 1.2 § f . 

In the next section, the results of these case studies will be presented and discussed. 


5. Discussions of Results 

We are interested in the effectiveness of the proposed global approach in minimizing the magnitudes of the 
base force and base moment as well as in comparing this approach with the local approach. We therefore, plotted 
the magnitudes of the base force and base moment (corresponding to the optimized trajectory) for each of three 
cases enumerated in the previous section. From an examination of the results shown in Figs. 6 through 11, the 
following important conclusions/points should be noted: 




1 - local 

2 - global 


Fig. 6 Base Reactions Profiles of Case Study (l),(a) Base Force Magnitude and (b) Base Moment 
Magnitude. 
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Fig. 1 1 Base Reactions Profiles for Case Study 4, (a) Base Force and (b) Base Moment. 
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(l)First consider the magnitude of the force transmitted to the base shown in Fig. 6(a) (for Case Study 1), Fig. 
7(a) (for Case Study 2), and Fig. 8(a) (for Gtse Study 3). Each one of these figures has two plots corresponding to 
the global and local approaches. As might be expected, the magnitude of the base force is the lowest for Case Study 
2 where only the base force is being minimized (see Fig. 7(a)). Furthermore, as shown in Fig. 7(a), the global 
approach "smooths" out the peak in the magnitude of the force which is seen in the local approach at t=0.8 second. 

4 (2) Next consider the magnitude of the moment transmitted to the base as shown in Fig. 6(b) (for Case Study 

1), Fig. 7(b) (for Case Study 2) and Fig. 8(b) (for Case Study 3). Once again each one of these figures has 2 plots 
corresponding to the local and global approaches. For the local approach, Case Studies (1) and (3) yield very small 
» base moment, while Case Study (2) (base force minimization only) yields a relatively high base moment. In 

contrast, for the global approach, the magnitudes of the moment transmitted to the base is virtually independent of 
the cost function (i.e. whether we are minimizing the magnitudes of the base force, the base moment or the 
combination of the two). We also see (Figs. 6(b), 7(b), and 8(b)) that the local approach is more effective than the 
global approach in minimizing the magnitudes of the base moment. 

(3) We can summarize the conclusions in (1) and (2) above by the following two statements: 

(a) The global approach is the better approach for reducing the magnitude 
of the base force whereas the local approach is superior for reducing the 
magnitude of the base moment. 

(b) In the global approach it is sufficient to minimize the magnitude of the 
base force. 

(4) The second derivative of the joint variables (joint "accelerations") for each one of the three joints are 
shown in Fig. 9, both for the local and global approaches: we observe that very similar acceleration profiles are 
obtained for these two approaches. The conepsonding manipulator configurations are shown in Fig. 10(a) (for the 
local approach) and Fig. 10(b) (for the global approach). One .should make careful note , both from Figs. 9 and 10, 
of the "sense" of motion of the three links. Title results shown in Figs. 9 and 10 correspond to Case Study 3. 

(5) The results of the simple sensitivity study are shown in Figs. 1 1(a) (magnitude of the base force) and Fig. 
1 1(b) (magnitude of the base moment) for the global approach and the cost function of Case Study 3. (The optimal 
results are shown in dashed lines with the perturbed results shown in solid lines). The following two conclusions can 
be reached from the results as shown in Fig. 11. 

1. The magnitude of the base moment is relatively insensitive to small changes in the trajectory 
(parameters). (Fig. 1 1 (b)). 

2. The magnitude of the base force is quite sensitive to 0 f , the final position of the redundant link (Fig. 

1 1(a)). (The "redundant" link was selected to be the link attached to the base in Fig. 10.) 

6. Summary and Conclusions 

In this paper we have developed a global approach for minimizing the base reactions of a kinematically 
redundant manipulator used in microgravity environments. The global approach utilizes the Partitioned Jacobian 
* method of redundancy resolution, the 4-3-4 method of joint-trajectory representation and the minimization of a cost 

function which is a time-integral of the magnitude of the base reactions. We also compared the global approach of 
this paper with a local approach developed earlier. 

The application of the global and local approaches to a planar 3 d.o.f. manipulator with 1 redundant degree of 
freedom demonstrated the effectiveness of the global approach in smoothing out undesired "peaks” in the force 
transmitted to the base which were observed in the local approach. We also learned that the global approach was 
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more effective in reducing the magnitude of the base force while the local approach was more effective in reducing 
the magnitude of the base moment. A simple sensitivity analysis performed on the optimal trajectory obtained for 
the global approach revealed that the magnitude of the force transmitted to the base is more sensitive to the final 
configuration of the links. We plan to do an extensive sensitivity analysis in the future. 

The approach will be implemented on a 4 d.o.f. traction drive manipulator for space application which is 
currently under construction. 
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Appendix 1 : Equations for 4-3-4 joint trajectory representation. 

The equations for 4-3-4 joint trajectory representation can be found in [10]. The equations of 0j(t), j=l,2,3 of 
the three-parameter representation are: 

For segment 1, 


0j (x) * (Of 0 a- a)x 4 + ox 3 + 0, 

(49) 

0 i(t)= [4(0 /o a '-Ox 3 + 3ax 2 ]/X! 

(50) 

e,(x) = [I2(e /tf a -a)x 2 +6ox]A 1 2 

(51) 


where 6 i = initial joint position, 

X = normalized time (0,l)=t/X l 

<*=f/g 

/=2e /o a'(6 + 2p'/p + p/pVe /o ap'/p(3 + 3'/p) + 2e /o a , p'/p 

g*pVp+3+3p/p'. 

For segment 2: 

02 (x )« (0 /o a - v , TyP-a, xf p 2 /2 ) x 3 + a, xf p^+VjXypx+Oy (52) 

0 2 (x) - [3(0 /o a , -v 1 x / p-aiX / 2 p 2 /2)x 2 +a 1 x / 2 p 2 x+v 1 x j P]A 2 (53) 

02(t) - [6(0 /o o'-v 1 x / p-a 1 x / 2 p 2 /2)x+a 1 x / 2 P 2 ]/x 2 2 (54) 

where 

x=0-x,)/x 2 

v,-0i(l)/x, 

a,«0,( l)/(P't/ 

For segment 3, 

03 (x) * (90 /o a'-4v 2 p'x / -a 2 x / 2 p' 2 /2)x 4 +(-80 /< ,a+3v 2 p / x / )x 3 + 
a 2 P 2 x / 2 x 2 /2+v 2 p'x / x+0 y 


( 55 ) 
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0 3 (x) = [4 (90 /o a -4v 2 p'y-o 2 x / 2 p <2 /:2)x 3 + 3(-80 /o a'+3 v 2 p' x f ) x 2 + 
a 2 p' 2 x / x+v 2 p x f ]/x 2 

(56) 

0 3 (x) = [ 12(90^,0 -4 v 2 p'x / r-a 2 x / 2 p' 2 ^)x 2 +6(-80 /o a'+3v 2 p' x f ) x + 
a 2 p' 2 x/]/r 3 2 

(57) 


where 

TsO-Tj-T^/Tj 

v 2 =30 / a/(x / P)-2v r a 1 Xy-P/2 
02=6 0 /o a/Cx^rp) 2 - 6 v, /(x 7 p ) - 2 o, 
0 2/ =0 2 (l). 


Appendix 2 : Base Reaction Dynamic Equations 

The following variables are defined 
m i = total mass of the link i 

/ ( = length of link i 

l ei = distance measured from the axis of joint (i-1) to the centroid of link i 
7j = moment of inertia about the centroid of link i 

For the 3 d.o.f. planar manipulator studied in this paper, the base force vector has a x-direction component, 
(F^ and a y-direction component, (F^y, and the base moment vector has only a z-direction component, (N 0 ) z . The 
base reactions dynamic equations are given by 

F b = M b q + V(q,q) (58) 

where F b =[(F o ) x (F 0 ) y (N o ) z ] T . 

Let the elements of M b be mb^, @=1,2,3) and (j= 1,2,3). And for convenience, we denote cos(i+...+j) by c t j 
and sin(i+...+j) by s L j. We have the following equations 


mZ, n =/ «3 w 3 J 123~ / 2 m 3' s 12 _/ e2 w 2' r 12 -, l( w 3 +m 2^r / cl w l 5 l 
mb\2~ ~l e i m 3 s \23~h m 3 s \2~^c2 m 2 s 12 
m *13® -/ c3 m 3 c 123 

mt >21= l c3 m 3 c 123 +l 2 m 3 c 12 +I c2 m 2 c 12 +l l( m 2 +m 3) c l +l cl m l c l 


mt>22 = ^ c 3 m 3 c 123 + ^2 m 3 c l2 + ^e2 m 2 c 12 


m *23 sl c2 m 3 c 123 


( 59 ) 
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m 6 3 1 - Vc3 m 3 c 23 + ^ Vc3 m 3 c 3 + 2 V2 m 3 + ^ Vc2 w 2 c 2 + *c3 2m 3 + 

/ 2 2 m 3 +/ j 2 «3+/ c2 2m 2 + ^ 1 2m 2 + ^c i i +^i 2 + ^3 

^32=^ / c3^ c 23 +2/ 2 / c3 m 3 c 3 +/ l / 2 m 3 c 2 +/ l / c2 m 2 c 2 +/ C 3 2m 3 +/ 2 2m 3 

+/ c2 2m 2 +/ 2 +/ 3 

m *33 ss ^l^c3 m 3 c 23 + ^2^c3 m 3 c 3 + ^c3 2m 3 c 3 + ^e3 2m 3 

and the elements of V(-) denoted by v ; are 

V,^[?3 2 +2(?2 + ?l)V?2 2+2 ?l?2 + ^l 2 ^c3 ,M 3 c 123 - t(?2 2+2 <?l<?2 + ^l 2 ) / 2 m 3- 

(?2 2+2 ?l?2 + ^l 2 ) / c2^ c 12 - <?l 2/ l m 3 + ?l 2/ l m 2 + ?l 2/ cl m l) c l 


v 2^[?3 2+ (^2 + 9l)?3 + ?2 2+ ^W2 + ?l 2 ^c3 m 3 5 123 - [^2 2+2 ^lWl 2 ) / 2 m 3 

-(?2 2 +2^ 1 ?2-?l 2 ) / c3^ J 12 - ^l 2/ l m 3 + ?l 2 ( / l m 2 +/ el m l^ 5 l 


V 3 =-[?3 2+2 ^2 + ^ 1 )?3 + ^2 2+2 <7 W2] 1 1 / c3 m 3 J 23'( 2 ^2^3 +2 ^1^3 + ^3 2 ) / 2 / c3 m 3- r 3 _ 

($2+24! ? 2 ) V2 m 3*3"(?2 +2 4 l^Vl / c2 m 2 i 2 C 60 ) 


Appendix 3:Physical Dimensions and Mass Properties 

The physical dimensions and the mass properties of the three degree-of-freedom planar manipulator are given 
in Table 1. The link length, mass, and inertia of the link are denoted by l^, and /, respectively. The centroid of 
the link is assumed to be at the center of the link. 


Table 1 : Simulation Data for 3 d.o.f. Manipulator 


link 

1 

2 

3 

m ( -(kg) 

1.0 

1.0 

1.0 

/,< m) 

0.5 

0.5 

0.5 

Ij(kg m 2 ) 

0.0208 

0.0208 

0.0208 


1 
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